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Abstract 



^ ■ The phase-ordering kinetics of the ferromagnetic two-dimensional Ising model with uniform bond 

c/3 . disorder is investigated by intensive Monte Carlo simulations. Simple ageing behaviour is observed 

in the single-time correlator and the two-time responses and correlators. The dynamical exponent 
z and the autocorrelation exponent Ac only depend on the ratio e/T, where e describes the width 



'^ . of the distribution of the disorder, whereas a more complicated behaviour is found for the non- 

a, 

O . equilibrium exponent a of the two-time response as well as for the autoresponse exponent A/j. The 

O 

scaling functions are observed to depend only on the dimensionless ratio e/T. If the length scales 

^ ■ are measured in terms of the time-dependent domain size L{t), the form of the scaling functions is 

^n : 

00 . in general independent of both e and T. Conditions limiting the validity of this 'superuniversality' 

o 
00 

o 



are discussed. 
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I. INTRODUCTION 

Aging phenomena have become one of the paradigmatic examples which are used to study 
fundamental aspects of non-equilibrium criticality, besides their practical use in materials 
processing. If some physical system is brought rapidly out of equilibrium by a sudden 
change of an external control parameter (an often-used device is to quench the system by 
lowering its temperature rapidly from the disordered high-temperature phase to its ordered 
low-temperature phase where at least two thermodynamically stable states co-exist) one 
often finds simultaneously the following three properties: (i) slow (i.e. non-exponential) 
dynamics, (ii) breaking of time-translation invariance and (iii) dynamical scaling, which 
together are said to constitute ageing behaviour. These features were first observed together 
in the mechanical properties of many polymeric materials by Struik [l| and it has since 
been understood that the broad characteristics of ageing can be found in many different 
types of non-equilibrium systems. Besides spin glasses |2|, y], other well-studied examples 
may be found in the phase-ordering kinetics of a ferromagnet quenched below its critical 
temperature T^ > ^ or in granular media |5|. 

The analysis of phase-ordering as it occurs in ferromagnets quenched from an initially 
disordered state into its coexistence phase with temperature T < Tc is particularly simple 
to formulate. The motion of the domain walls is driven by the surface tension between the 
ordered domains. The typical time- dependent length scale is related to the linear size of 
these ordered domains which grows as L = L(t) ~ t^/^ where z is the dynamical exponent. 
For a non-conserved order-parameter, it is well-known that 2 = 2 [j] . Because of the simple 
algebraic scaling of the linear domain size L{t), one expects the following scaling behaviour 
for the single-time correlation function 



C{t; r) := (0(t, r)0(t, 0)) = r'F (^-^J (l) 

for sufficiently large times t ^ tmicro (with a microscopic reference time tmicro such that 
-^(^micro) IS of the ordcr of the lattice constant). Similarly, for the two-time correlation and 
response functions (in the ageing regime, where the observation time t and the waiting time 
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Here 0(t, r) is the space-tiine-dependent order-parameter, whereas h{s, r) is the conjugate 
magnetic field (spatial translation-invariance of all averages will be assumed throughout this 
paper) and a and b are ageing exponents. The scaling functions Fcfi{y, 0) ~ y-^cR/^ for 
y ^ oo which defines the autocorrelation exponent Xc and the autoresponse exponent Xr. 
For phase-ordering kinetics in pure systems, it is generally admitted that 6 = and simple 
scaling arguments show that a = 1/z. For an initial high-temperature state and for pure 
ferromagnets, Ac = A/? is independent of the known equilibrium exponents [2|, U, [y, |7|, |8|, |9|, 



lOj. The scenario just described is referred to as simple ageing. The conditions for the onset 



of ageing and the possible scaling forms have been carefully discussed in UJ, [l2| . 

In going from phase-ordering kinetics in simple ferromagnets to glassy systems (usually 
modelled by spin systems with disorder and frustration), one expects more complicated 
growth laws L = L{t) which describe a cross-over from the domain growth of essentially pure 
systems (as long as L{t) is small compared to the typical distance between disorder-created 
defects) to a late-time regime with a slower growth and dominated by the defect structure. 
For disordered ferromagnets without frustration, this can be studied through generalisations 
of the AUen-Cahn equation, which attempt to describe how the pinning of the domain walls 
created by the disorder should be overcome by thermal activation [13|, \l^ . In addition, it 
was suggested that once the unique reference length scale is chosen to be L{t), the resulting 
scaling functions should become superuniversal in the sense that they should be independent 
of the disorder 15|. This sup eruni versality has indeed been confirmed for the single-time 



spin-spin correlator 
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2l|. Furthermore, superuniversality is in qualitative 



agreement with the experimental observation in several distinct polymers and metals that 
the linear response to a small mechanical stress can be described in terms of an universal 
master curve which is independent of the material studied p| . The universal sca ling functions 
of the pure case can be calculated from the theory of local scale-invariance [2^. Tests in 



non-integrable systems inc 



ude the auto- and space-time-responses in the two- and three- 



dimensional Ising models 23|], the autocorrelation function in the two-dimensional Ising 



model [2J|, the same quantities in the two-dimensional g-states Potts model with g = 2, 3, 8 
25l | as well as the autoresponse function in several other cases, see 26| for a recent review. 
Therefore, with the help of superuniversality, if confirmed, the scaling functions describing 
non-equilibrium relaxation of quite complex systems would become analytically treatable. 

In this work, we shall consider a two-dimensional ferromagnetic Isin g m odel with 
quenched bond disorder. The nearest-neighbour Hamiltonian is given by 271. |28| 
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The random variables Jij are uniformly distributed over [1 — e/2, 1 + 6/2] where < e < 2. 
The model has a second-order phase transition at a critical temperature Tc{e) > between 
a paramagnetic and a ferromagnetic state. It is thought that Tc{e) ~ ^c(O) — 2.269 . . . 
should not depend strongly on e. Using heat-bath dynamics with a non-conserved order- 
parameter and starting from a fully disordered initial state, phase-ordering occurs and there 
is evidence which suggests that the characteristic length scale L = L(t) ~ t^'^ should scale 
algebraically and where the dynamical exponent z = z{T, e) should depend continuously on 
;he temperature and the disorder e. Indeed, generalising the Huse- Henley heuristic argument 
13| by considering the case when the disorder-created energy barriers for the motion of the 
domain walls are distributed logarithmically with respect to the domain size L{t), the form 

z = z{T, e) = 2 + e/T . (5) 

was proposed, where the constant e parametrises the barrier height. Simulations of the linear 
domain size 27| seemed to confirm this, with the empirical identification e = e and are 



also consistent with the results of field-theoretical studies in the Cardy-Ostlund model 30|. 
Data from the thermoremanent magnetisation Mtrm(^, ■§) = h j^duR(t,u) = s~°'fM(t/s) 
were used to estimate the exponent a and, assuming a = 1/z, also looked consistent with 

n 

([5]) and e = e, at least for values for T and e for which z did not become too large [31| 
Superuniversality has been confirmed recently for the 
hand, the conclusion of a simple ageing reached in [ 
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more recent simulations for the random-site Ising model 33] . In that work, a scaling form 
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321 ]. On the other 



291 has been questioned by 



C{t, s) = C,t{t -s) + C^g{h{t - s)/h{s)) with h{u) 



exp 






was considered, where /i is 



a fit parameter. In the limit /i — ^ 1, one recovers the simple ageing scenario described above 



(the stationary part Cst{t—s) merely represents an irrelevant correction to the leading scaling 
behaviour). The case /i < 1 is called sub-ageing, and the case /i > 1 is called super- ageing. 



In (33|, systematic deviations from the dynamical scaling of simple ageing were observed in 
the 2D random-site Ising model quenched to below T^. A data collapse could be achieved, 
however, by allowing /i to vary and values in the range /i ~ 1.03 — 1.04 were reported 33|. 



This finding was interpreted as to suggest the presence of a slight super-ageing effect (33j. 

In the following, we shall present new data on the single-time and two-time correlations 
as well as on two-time response functions. As we shall see in section II, our data are fully 
compatible with the simple ageing scenario and furthermore, looking at a larger range of 
values of z, we find that the dynamical exponent z = z{T,e) = z{e/T) depends on the 
control parameters in a more complicated way than suggested in eq. ([5]). These conclusions 
are also valid for the two-time response function. The dependence of the various non- 
equilibrium exponents on both e and T will be studied. We also show evidence that the 
scaling functions only depend on the ratio e/T of the control parameters and finally confirm 
the generic superuniversality of the scaling functions of correlation and response functions. 
However, we also find two conditions which must be satisfied for superuniversality to hold. 
Our conclusions are given in section III. 

II. RESULTS 

The simulations are carried out as follows. For the integrated response we simulated 
systems with A/" = 300^ spins using the standard heat-bath algorithm. Prepared in an 
uncorrelated initial state corresponding to infinite temperatures, the system is quenched to 
the final temperature in the presence of a random binary field hi = ±h with strength h = 



0.05, following the well-established method of Barrat [34] (using a random field avoids a bias 
which would drive the system rapidly out of the scaling regime). Turning off the random field 
after a waiting time s, the thermoremanent magnetisation Mtrm(^5 s) = t^ Yli {^ii't)hi{s)) 
is measured at time t. We averaged over at least 5-10^ different runs with different initial 
states and different realizations of the noise. We point out that the data discussed in this 



paper are of much higher quality than our earlier data 3l| for the autoresponse. For the 
autocorrelation function, we considered systems with up to 600^ spins in order to avoid the 
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FIG. 1: Scaling of the autocorrelation C{t,s) for (a) e = 0.5, T = 1, (b) e = 1, T = 0.4 and (c) 
e = 2, T = 0.4, for several values of s. Here and in the following error bars are smaller than the 
sizes of the symbols, unless explicitly stated otherwise. 

appearance of finite-size effects for the times accessed in the simulations. The data discussed 
in the following have been obtained after averaging over at least 5000 different runs with 
different random numbers. Our main focus was on e = 0.5, e = 1, and e = 2 where we 
considered for every case at least four different temperatures. In addition, some runs where 
also done for other values of e. The total study took approximately 2 ■ 10^ CPU hours 
on Virginia Tech's System X supercomputer composed of Dual 2.3 GHz PowerPC 970FX 
processors. 

Our first question is about the scaling form to be used. In figure [1], we show data 
for the autocorrelation C{t,s) = A/"^^ ^,-(crj(t)a"j(s)), plotted over against t/s, for several 
typical values of e and T. A nice data collapse is seen, which is fully consistent with simple 
ageing. Our scaling plots also imply that the exponent 6 = 0, analogously to what is 
found in the phase-ordering of pure systems. In 33|, a 'super-ageing' scaling form C(t, s) = 
C(exp[ ^ ~^' ^_ ~^ — ]) was considered for the random-site Ising model where the exponent 
/i > 1 is fitted to the data. Simple ageing is recovered in the // — *> 1 limit. However, the 



values of /i ~ 1.03 reported in 



are so close to unity that a careful study on possible 
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FIG. 2: Left panel: Dynamical exponent z, determined from the condition C{t,L[t)) = 21 as a 
function of e/T. The dashed line corresponds to the equation ^. Right panel: ageing exponent 
o, as determined from the scaling of MtrmI^, s). In some cases we have more than one data point 
for a given value of e/T. These data points correspond to different values of e and T with e/T 
constant, see table HI showing that a is not simply a function of the ratio e/T but that it depends 
in a more complicated way on both e and T. 

finite-time corrections to scaling appears to be required before such a conclusion could be 
accepted 35||. For all e < 2, our data for the random-bond model show no hint for a 'super- 



ageing' behaviour [37[, in contrast to the findings in [i 



Note that the observed scaling 



form of simple ageing would be incompatible with a non-power-law form of L{t) in the range 
of times considered. We also remark that for e = 2 and larger temperatures finite-time 
corrections to simple ageing are observed, see 



31 



36|. 



Having in this way checked that the relevant length scale L = L{t) should indeed scale 
algebraically with time, we next determined the dynamical exponent z = z(T,e) from the 



criterion 



involving the single-time correlator 

1 



C(t,L(t)) 



, L(t)~t 



l/z{T,e) 



(6) 



The results for z are shown in figure [2] and listed in table [H First, we observe that the values 
of z obtained in fact only depend on the dimensionless ratio e/T, to within our numerical 



e T 
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a 
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Ai?/^ 


0.5 1.5 


2.08(1) 
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0.59(1) 
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0.5 1.0 


2.11(1) 


0.41(1) 


0.565(10) 


0.61(2) 


0.5 0.8 


2.16(1) 


0.41(1) 


0.56(1) 


0.61(3) 


0.5 0.6 


2.33(1) 


0.40(1) 


0.54(1) 


0.60(2) 


0.5 0.5 


2.46(1) 


0.36(1) 


0.48(1) 


0.55(2) 


0.5 0.4 


2.64(2) 


0.33(1) 


0.46(1) 


0.50(2) 


0.5 0.3 


3.02(2) 


0.29(1) 


0.385(10) 


0.46(2) 


1.0 1.0 


2.45(1) 


0.32(1) 


0.49(1) 


0.51(2) 


1.0 0.8 


2.65(2) 


0.31(1) 


0.445(10) 


0.50(2) 


1.0 0.6 


3.02(2) 


0.25(1) 


0.38(1) 


0.43(2) 


1.0 0.4 


3.85(3) 


0.17(1) 


0.29(1) 


0.34(1) 


1.5 0.9 


3.02(3) 


- 


0.375(10) 


- 


2.0 1.0 


3.39(3) 


0.10(1) 


0.315(10) 


0.35(2) 


2.0 0.8 


3.92(4) 


0.09(1) 


0.270(5) 


0.32(2) 


2.0 0.6 


4.97(4) 


0.075(10) 


0.217(3) 


0.28(1) 


2.0 0.5 


5.76(5) 


- 


0.189(3) 


- 


2.0 0.4 


7.22(6) 


0.05(1) 


0.155(3) 


0.21(1) 



TABLE I: Dynamical exponent z, non-equilibrium response exponent a, the autocorrelation expo- 
nent Xc/ z and the autoresponse exponent Xr/ z for different values of e and T. 



accuracy. Second, we see that the function z = z{e/T) is non-hnear and only becomes an 
approximately linear function in a relatively small region of values of z. We are confident 
that our results are more reliable than earlier ones since they do not just describe the scaling 
of a single quantity, but rather will be needed for the correct scaling description of several 
other observables, as we shall show below. We stress that only the values of z as given in 
table [I will lead to a good scaling according to simple ageing without having to consider 
possible corrections to scaling. 

In the same way, in figure [3] we show data for the scaling of the thermoremanent mag- 
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FIG. 3: Scaling of the thermoremanent magnetisation MtrmI^j s) for several values of s, with (a) 
e = 0.5, r = 0.4, (b) e = 1, T = 0.6 and (c) e = 2, T = 0.8. 



netisation, expected to be of the form 



MTRM(t,s)=S-VM(t/s) . 



(7) 



The exponent a is obtained in the usual way 39| by plotting the thermoremanent magneti- 
sation as a function of the waiting time s for fixed values of the ratio tj s. The resulting 
power-law decay, see eq. (jTj), then yields the value of the exponent a. The numerical values 
and error bars given in table [T] and shown in figure [2] are obtained after averaging over the 
values obtained for five different values of t/s, namely 3, 5, 7, 10, and 15. Looking at these 
values (see figure [2]), we observe that the estimates a = a{T,e) scatter considerably more 
than those for z. We consider this scatter to be large enough to conclude that a cannot be 
reduced to a function of the single variable e/T. Furthermore, considering in detail the nu- 



merical values from tab 
of pure ferromagnets 



elU we see that the relation a = 1/z, known from the phase-ordering 



4l|, is no longer valid. 



In order to understand this finding, let us briefly reconsider how the relation a = 1/z may 
be derived for pure ferromagnets. Consider a pure ferromagnet in an external oscillating 
magnetic field of angular frequency u. The dissipative part of the linear response is given 



as the imaginary part of the dynamic susceptibihty and reads (see e.g. |2|) 

X"{u;,t)= /dni?(t,n)sin(^(t-n)) =XiM + rV2M) + --- (8) 



where the last relation follows from the usually assumed scaling ([3]) of the autoresponse 
function R{t, s). On the other hand, motivated from the physical picture that the dynamics 



rom the motion of the domain walls between the ordered 



in phase-ordering should only come 
domains, one would expect to find 40| 

x"(a;,t) = x:tM + ^W-V4eM) + ... (9) 

from which one may identify the stationary and the ageing part with the terms in eq. ([8]) 
coming from the scaling analysis. Since only the domain boundaries contribute to the 
dynamics, the leading time-dependent part should be proportional to the surface area of 
the domain divided by the total volume, hence to L{tY~^ / L{tY = ^/ L{t) which accounts 
for the factor 1/ L{t) in Q. Comparison of eqs. (151191) . together with L{t) ~ t^/^, then gives 
a = 1/z. 

Our empirical observation that az < 1 suggests that the above argument should no longer 
apply to random ferromagnets. Since ([H]) only depends on the dynamical scaling assumption 
([3]), and given that our numerical results appear to be compatible with it, we expect that 
(IH]) should remain valid for disordered ferromagnets. Since also for disordered ferromagnets, 
the contribution to the ageing behaviour should come from the boundary region between 
ordered domains (this is also suggested by looking at the microscopic spin configurations. 



see e.g. (36|]), we expect it to be proportional to Nd{L)/Nb{L) where Nb^d{L) denote the 
number of mobile spins in the bulk and at the domain boundaries, respectively. While one 
should still have Nb{L) ~ L'^, disorder may cause the domain boundary to become fractal 
and, hence, Nd{L) ~ L*^* with df the fractal dimension (and d{ = d — 1 for the pure case). 
Then eq. Q) would be replaced by 

X"(u;, t) = x'Li^) + m-^'-'^^^x'UuJt) + . . . (10) 

and comparison with ([8]) would now imply 

d-di 
a = ^. (11) 

Our empirical results (table [T]) imply that df > d — 1, that is, the disorder should modify 
the domain boundaries into fractal curves. From eq. (ITTl) . since a depends on the dynamical 

10 



1.0 



b 

*. o t=500 1 
Q ■ t=1000 
o t=2000 




0.5 1 1.5 2 
r/L(t) 




' 1 


1 ' 1 ' 1 ' 




o s=100 




■ s=200 ■ 




* s=400 


"(c)* 

1 1 


■ 



0.5 1 1.5 2 
r/L(s) 



12 3 4 5 
r/L(s) 



FIG. 4: (a) Scaling behaviour of the single-time correlator C{t, r) for e = 1 and T = 0.8 and several 
values of t, as a function of the scaling variable r/L{t). (b) Scaling of the space-time-dependent 
correlator C(t, s; r) for e = 0.5 and T = 0.5 for several values of s and two fixed values of y = t/s, as 
a function of the scaling variable r/L{s). (c) Scaling of the space-time-dependent thermorenianent 
magnetisation MTRM{t-, s; r) for e = 2 and T = 0.6 for several values of s, with y = t/s = 5. The 
lines are guides to the eye. 

exponent z = z{e/T) as well as the fractal dimension d{, it may appear more natural that a 
cannot be written as a function of the single variable e/T. 

In figure H] we show the scaling behaviour of the space- and time-dependent correlation 
and response functions for various values of the waiting time s. For the selected typical 
values of e and T, a simple ageing behaviour is observed, in agreement with the observed 
scaling behaviour of the autocorrelation and of the autoresponse. We have found completely 
analogous results for all other values of e and of T which we considered. In fact, the 
scaling is much cleaner than for the r = quantities and for none of the studied cases 
a sizeable correction to scaling could be identified. Obviously, space- and time-dependent 
quantities are much better suited for an investigation of the scaling forms than quantities 
that only depend on time. Similar conclusions have recently been drawn from a study of 
nonequilibrium growth models 421 ]. 
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FIG. 5: (a) Space-time correlator C(t, s; r) for e = T = 1 (filled symbols) and e = T = 0.5 (open 
symbols) for several waiting times s and the values y = t/s = 2 and y = t/s = 6. The full curves 
are guides to the eye. (b) Autocorrelator C{t, s) for e = T = 1 (filled symbols) and e = T = 0.5 
(open symbols) for several waiting times s. (c) Scaled space-time thermoremanent magnetisation 
for e = T = 1 (filled symbols) (with a = 0.32) and e = T = 0.5 (open symbols) (with a = 0.36) for 
several waiting times s and the values y = t/s = 2 and y = t/s = 6. The full curves are guides to 
the eye. 

Using the scaling forms (12|3|) in the limit of large y = t/s for the autocorrelation and 
autoresponse functions (where r = 0), we have also extracted the exponents \c/ z and \r/ z 
and list our results in table [B In contrast to the pure case, where for fully disordered initial 
conditions one may show that Ac = \r |4, 11Q| , the values of the autocorrelation exponent 
\c are different from those of the autoresponse exponent \r. In particular, we find that 
within our numerical accuracy, \c/ z is a function of the single variable e/T, at least for 
e < 2, while \r/ z cannot be expressed in this way. Our data suggest that Ap > A^ and they 
are consistent with the rigorous Yeung-Rao-Desai inequality Ac > d/2 43|. Furthermore, 
we observe that \r/ z — a should be practically constant (again for e < 2 and with a value 
in the range ~ 0.17 — 0.20). 

In conclusion, our data are clearly consistent with simple ageing of the single-time and 
two-time correlation functions, as well as for the thermoremanent magnetisation, and fully 
confirm the anticipated scaling forms (TT][2|3|) with L(t) ~ t^/^ . 
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Next, we shall compare the form of the scaling functions, for several values of the control 
parameters e and T. Since the dynamical exponent z = z{e/T), see figure [21 one might 
expect that the scaling functions themselves should only depend on the ratio e/T, rather 
than on e and T separately. In figure we test this idea by comparing data for e = T = 1 
with those for e = T = 0.5, for the three cases of (a) the space-time-dependent correlation 
C(t, s; r), (b) the autocorrelation C{t, s) and (c) the space-time-dependent thermoremanent 
magnetisation Mtrm(^, ■5;r) = h J^duR(t,u; r). In all cases, there is a clear scaling be- 



haviour consistent with simple ageing [4J] and the scaling functions nicely superpose (for 
the integrated response, the data for e = T = 0.5 were multiplied by 1.17 in order to take 
into account the well-known presence of the non-universal numerical prefactor). This result, 
namely that the form of the scaling function only depends on the ratio e/T, goes beyond 
the standard scaling form flT][2][3l) . yet it does not require to re-scale the length by the typical 
domain size L{t), as it would be required for a test of superuniversality. More systematic 
tests of this result would be welcome. We point out that the findings of figure [5] are consis- 
tent with our earlier observation that Ac = \c{s/T). The more complicated dependence of 
\r on both e and T would only appear if in plots such as figure [5b one would concentrate 
on the region \r\ k, 0. 

After these preparations, we are ready for a test of superuniversality. Superuniversality 
of the single-time correlator C(t, r) is tested in figure [6l First, we show in figure [6^ data for 
several values of < e < 2 and T, where the values of z may be read off from figure [2] and 
table [T| The times were chosen sufficiently large such that a clear scaling behaviour has set 
in. Using the typical length scale L{t) as determined earlier from eq. ([6]), we see that indeed 
all data, including the one for the case without disorder, collapse nicely onto a single curve. 



within the numerical accuracy. This is a a 



m y, m 



ear confirmation of superuniversality, very much 



in agreement with earlier studies 16|, [l7|, [l8|, [l9|, [20|, |2l|, [3^. However, when considering 
the case £ = 2, which is shown in figure [6)d, a different picture emerges. Clearly, the scaling 
curves for C(t, r) as obtained for several values of T again collapse onto each other, but, as 
the comparison with the scaling function of the pure case shows, the scaling function is no 
longer the same as the one of the pure case. Therefore, if taken at face value, the case e = 2 
might represent a distinct superuniversality class. In order to get a better understanding 
on this point, we show in figure [Gb data for a relatively large value of e. One sees that 
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E=0.5,T=1.0, t=9600 
E=0.5, T=0.4, t=600 
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E=0, T=1.0, t=600 . 




1.0 
0.8 
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E=2, T=0.8, t=2400 
E=2, T=0.6, t=4800 
E=2, T=0.4, t=9600 
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t=60 
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t= 12000 



r/L(t) 



r/L(t) 



r/L(t) 



FIG. 6: (Colour online) Test of the superuniversality for the single-time correlation function. In 
panel (a), data for e = 0,0.5, 1, 1.5 are shown to fall onto the same rescaled curve. In panel (b), 
we show that the data for e = 2 describe a curve different from the one of panel (a). The green 
curve is the scaling function of the pure system. In panel (c), data for e = 1.6 and different times 
are compared with the scaling curves as measured for e = (lower green line) and e = 2 (upper 
green line). 

for moderately large times, the data are very close to the curve found for e = 2, but when 
the time t is made very large, a cross-over to the curve of the pure case e = is observed. 
Qualitatively, the cross-over time tx increases when e — »■ 2 and becomes so large that a cross- 
over is no longer detectable for the times (and the lattice sizes, which must be increased to 
large times in order to avoid finite-size effects) reachable with our numerical methods. We 
did not see any sign for a cross-over in our data with £ = 2, but purely numerical techniques 
cannot distinguish between a very large and an infinite cross-over time tx- 

We extend the test of superuniversality to the case of space-time-dependent two-time cor- 
relators in figure [7] and similarly for the space-time-dependent two-time response in figure O 
Qualitatively, we arrive at essentially the same conclusion as for the single-time correlator. 
In both cases, panel (a) demonstrates a superuniversal behaviour for values of < e < 2. 
We show here our data for y = t/s = 4, but the same behaviour is observed for other values 
of t/s accessed in this study, namely 2 < t/s < 10. However, closer inspection also shows 
that superuniversality is no longer true for relatively small spatial distances \r\/L{s) < 0.5. 
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• 8=2.0, T=1.0 
□ £=2.0, T=0.8 

♦ £=2.0, T=0.6 
O £=2.0, T=0.4 
£=0,T=1.0 




FIG. 7: (Colour online) Test of the superuniversality of the two-time correlation function, with 
y = t/s = 4 and analogous to Figure [H The data shown are for s = 100. The green curves show 
the data of the pure system. The inset in panel (a) is a magnification of the area within the dashed 
box. The black lines are guides to the eye. 

This is shown in the insets of the panels (a) in figures [7] and [H Indeed, systematic devia- 
tions are observed for small spatial distances, the largest deviations being observed for the 
autocorrelation and autoresponse functions with \r\ = 0. The value of \r\/L{s) where the 
deviations set in seems to depend slightly on the value of s, but a larger range of s values 
than accessed in the present study is needed for a more quantitative discussion of this point. 
That means that although dynamical scaling does hold true even down to the autocorrela- 
tors and autoresponses, correlators and responses taken over a spatial distance of at least a 
typical cluster size L{t) show yet a larger degree of universality. 

This requirement appears to be consistent with the known numerical values of the expo- 
nents Xc/z and Xr/z as listed in table [B Superuniversality at r = would have required 
that their values should have been equal to those of the pure case e = 0, but we have rather 
seen that they depend on e and T. 

On the other hand, the case e = 2 again stands apart, as we illustrate in the panels (b) 
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• E=2.0, T=1.0 
□ E=2.0, T=0.8 

♦ E=2.0, T=0.6 
£=0, T=1.0 




FIG. 8: (Colour online) Test of the superuniversality of the two-time thermoremanent magnetisa- 
tion, with y = t/s = 4 and analogous to Figure [6j The data shown are for s = 100. In order to 
compare the data for different e and T, we have to take care of the fact that the response contains 
a non-universal multiplicative factor. The data have been multiplied by a constant in such a way 
that the different curves coincide with the curve for e = and T = 1 for large values of r/L{s). 
The green curves show the data of the pure system. The inset in panel (a) is a magnification of 
the area within the dashed box. The black lines are guides to the eye. 

in both figure [7] and [HI Comparing the data for e = 2 with the scaling functions found for 
< 5 < 2 (see panels (a)), we find small but systematic deviations. We stress that although 
these deviations are not very large, they are well outside the error bars of our high-quality 
data. As for the single-time correlator, it remains a possibility that for enormous times there 
might occur a cross-over to the scaling function of the pure case, but the relevant cross-over 
time tx is far larger than the time scales reached by our simulation. 
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III. DISCUSSION 

In this work, we have studied the non-equihbriuin scahng behaviour of a disordered 
Ising model without frustration, in an attempt to appreciate better the role of disorder by 
considering its effects in a system which is no longer identical to a pure system but which 
yet does not show the full complexity of a spin glass. Our conclusions are as follows: 

1. When quenched to a temperature T < T^ from a totally disordered state, the two- 
dimensional bond-disordered Ising model undergoes phase-ordering kinetics. The typ- 
ical length scale L{t) ~ t^/^ of the ordered domains scales algebraically with time, 
where the dynamical exponent z = z{e/T) depends continuously on the dimensionless 
ratio e/T of the control parameters. 

Quantitatively, this dependence can be read off from figure [2], and does not agree with 
earlier proposals of a linear relation such as in eq. ([5]). We have seen that this value 
of z correctly describes the dynamical scaling of not only single-time correlators, but 
of the two-time correlators and responses as well. 

2. Our data are completely compatible with simple ageing. 

3. While the non-equilibrium exponent 6 = of the correlation function is un- modified 
with respect to pure phase-ordering systems, the non-equilibrium exponent a which 
describes the scaling of the response function is no longer simply related to the dy- 
namical exponent z, see table [T] and figure El 

We propose to account for this finding in terms of a postulated fractal structure of 
the domain walls, which has led us to eq. (TTTj) . Further tests of this idea would be 
welcome. 

The autocorrelation and autoresponse exponents Ac = \c{s/T) and \r = Xfi{T,e) 
are distinct from each other, in contrast to the pure case. 

4. Our data suggest that the form of the scaling functions only depends on the dimension- 
less ratio e/T. It remains to be seen to what extent this observation can be extended 
to different systems. 
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5. In general, our data appear to confirm the superuniversality hypothesis, that is when 
all length scales are expressed in terms of L{t), the form of the scaling function is 
independent of both the disorder e and the temperature T. 

However, we have also found two important qualifications: 

(a) Superuniversality does not hold for sufficiently small spatial distances |r'|/L(t) < 
0.5. 

(b) For e = 2, although we find throughout a similar data collapse, the form of the 
scaling functions no longer coincide with the ones of the pure case. It is not 
understood whether the data presented here should be viewed as giving evidence 
for a distinct superuniversality class or else if there is a cross-over to the scaling 
functions of the pure scale at time scales much larger than the ones reached in 
our study. 

A better understanding of superuniversality will require an explanation of these con- 
ditions. 

What can these findings tell us on the behaviour of real materials ? Indeed, it has 
been shown recently, in a comparative study of the three-dimensional lajidoni-field Ising 



model and the three-dimensional Edwards- Anderson spin glass 2l|, that superuniversality 



is apparently satisfied in the former case (in which the disorder is 'weak' such that the ground 
state is still ferromagnetically ordered) while in the latter it is not 45|. Our own result is 
in qualitative agreement with this, but it raises the question how to explain the celebrated 
universality of the scaling functions for the linear response found in largely different materials 
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